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Abstract 

Quantum Brownian oscillator model (QBM), in the Fock-space representation, can 
be viewed as a multi-level spin-boson model. At sufficiently low temperature, the 
oscillator degrees of freedom are dynamically reduced to the lowest two levels and 
the system behaves effectively as a two-level (E2L) spin-boson model (SBM) in this 
limit. We discuss the physical mechanism of level reduction and analyze the behavior 
of E2L-SBM from the QBM solutions. The availability of close solutions for the QBM 
enables us to study the non-Markovian features of decoherence and leakage in a SBM in 
the non-perturbative regime (e.g. without invoking the Born approximation) in better 
details than before. Our result captures very well the characteristic non-Markovian 
short time low temperature behavior common in many models. 
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1 Introduction 



Recent development in quantum information processing and quantum computation has at- 
tracted much attention to the study of discrete quantum systems with finite degrees of 
freedom. The most commonly used model is an array of interacting two-level systems (2LS) 
each of which representing a qubit. As the system almost always interacts with its envi- 
ronment, quantum decoherence in the system usually is the most serious obstacle of actual 
implementation of quantum information processors El IS] ■ For this reason a detailed un- 
derstanding of quantum decoherence in open systems is crucial. There are a handful models 
useful for such studies, the quantum Brownian motion(QBM) jH E] is one, the spin- 
boson model jH [H] is another: the system in the former case is a harmonic oscillator and in 
the latter case a 2LS, both interacting with an environment of a harmonic oscillator bath 
(HOB). 

Most qubit models presently employed are the results of picking out the levels most 
relevant to the description of the qubit from a multi-level structure. In atom optics, internal 
electronic excitations are often approximated by a 2LS consisting of the ground state and 
the excited state. A similar model is used for the study of low temperature tunneling 
process where the two levels degrees of freedom represent the quasi ground states in a double 
well potential. The simplification to two-levels allows for detailed analytical or numerical 
treatment, but this remains an approximation applicable only when the effect of higher 
levels are negligible, e.g. , at low enough temperature when higher levels are not populated. 
However, in the presence of gate operation, the existence of higher levels causes a leakage 
of the 2LS due to transitions to other levels. Some extra perturbation may be necessary 
to select or restrict multi-level structure into the particular levels of interest [H|. In order to 
make a quantitative estimation of decoherence with a leakage effect, it is more desirable to 
study open system models which maintain the multi-level structure. 

In the present paper, we study the aspects of realistic qubits naturally arising from QBM, 
taking advantage of a fairly good understanding from the detailed studies over the last few 
decades. In particular, we focus on harmonic QBM, which can be viewed as an oo-level spin- 
boson model. Commonly used two-level spin-boson model can be obtained by restricting 
the harmonic oscillator Fock space to the lowest two levels. This correspondence allows for 
a detailed analysis of the spin-boson model from the known results of QBM. In particular, 
we will focus on the non-Markovian aspects of decoherence. Non-Markovian dynamics, 
often neglected in the literature (models are mainly based on a Markov approximation) for 
technical simplicity, is actually of crucial importance for the realistic implementations of 
quantum information processing. The 'effective' model we consider here invokes a two level 
simplification from a multi-level structure. How realistic this is certainly depends on the 
way the qubits are defined and realized in the multi-level structure usually encountered in 
actual experimental conditions. Nevertheless our model is able to capture the characteristic 
short time behavior common in many physical examples. 

Beyond the commonly assumed Ohmic spectrum for the bath, generic non-Ohmic en- 
vironments can be studied with this model. Contrary to the Ohmic case, the sub-Ohmic 
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environment (including l/f a type) causes nontrivial long time behavior such as anomalous 
diffusion or localization [H] owing to the long range temporal bath correlation. As demon- 
strated in ^0], the influence of slow environment can be dynamically decoupled from the 
system by using relatively slow pulses [TT]. In the present paper, we will mainly focus on the 
opposite case of supra-Ohmic environments [U El EE] Owing to the ultra-short time bath 
correlations, nontrivial short-time system dynamics enters, which is particularly difficult to 
describe by means of other models or approximations. The decoherece time scale in the 
supra-Ohmic environment can be much shorter than the one in the Ohmic case and thus is 
hard to remove by external pulses. Thus supra-Ohmic environment can be a major obstacle 
for the realization of quantum computation and information processing. 

We need to emphasize that to fully follow the coherence of an open system where self- 
consistency is required, because of the back-action from the environment, we need to study 
non-Markovian processs. We will also argue that, for a generic class of environment, Marko- 
vian approximations are not strictly valid. To facilitate comparison with results in related 
papers we will compare our methods with other commonly used approximations to the 
spin-boson model, such as the Born approximation and the Born-Markov rotating-wave 
approximation^] for two-level and multi- level systems. 

The outline of this paper is as follows: In Sec. El we specify the model and cast it in the 
influence functional formalism in the presence of an external field. In Section El we outline 
our idea of an effective 2LS using the QBM approach. Then we make correspondence be- 
tween the phase space representation discussed in Sec. El with the Fock space representation. 
We compare our approach with other methods based on Born-Markov and rotating-wave- 
approximation. Our results are presented in Sec. 14.11 In Sec. 14.21 we discuss the limitations 
and potential extensions of this approach. 

2 QBM in the presence of an external field 
2.1 The model 

Our model consists of a Brownian particle interacting with a thermal bath in the presence 
of an external field. We follow the notion developed in EI]. (We use the units in which 
ks = h = 1.) The Hamiltonian for this model can be written as 



where the dynamics of the system S (with coordinate x and momentum p) is described by 
the Hamiltonian 



and the (bare) potential V (x) is related to the physical potential by a counter term AV^ i.e. 
V (x) — V(x) + AV (see below). The Hamiltonian of the bath is assumed to be composed 



H = Hs + Hb + Hj + Hp, 





(2) 



3 



of harmonic oscillators with natural frequencies oo n and masses m n , 

N 2 2 2 

H b = 1^(^ + g } - (3) 

n=l Z 

where (gi, qN,Pi, ■■■,Pn) are the coordinates and their conjugate momenta. The interac- 
tion between the system S and the bath B is assumed to be bilinear, 

N 

Hl = X^ CnQn, (4) 
n=l 

where c n is the coupling constant between the Brownian oscillator and the nth bath oscillator 
with coordinate q n . The coupling constants are related to the spectral density Jb(cu) of the 
bath by, 

J fl ( w ) = 7r£-^— (5) 

We asssume the spectral density has the form 

J b (lo) = 2M 1 uj u e- u/A , (6) 

where v — 1 is Ohmic, v < 1 is sub-Ohmic , and z/ > 1 is supra-Ohmic. We will discuss the 
Ohmic and supra-Ohmic [y = 3) cases in detail. 

The counter term AV depends on c n , m n , u n , p and x and is given by 

J 2M 7 Ax 2 /7r (z/=l) 

\ 2M 1 Ap 2 /M 2 ir + 2M 7 AV/7T (^ = 3). l J 

This term is introduced to cancel the shift in the mass and frequency of the Brownian 
oscillator due to its interaction with the bath which will become divergent when the frequency 
cutoff A — > oo. As is customary, we consider the renormalized quantities after including a 
counter term as the physical observables with specified values. 
For a linear QBM, the potential is 

j r / s mvi 2 x 2 

V{x) = —j— , (8) 

where Q is the natural frequency of the system oscillator 
Finally, the Hamiltonian for the external field is 

Hp = -xE(t), (9) 

where E(t) is the external field. 
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2.2 The Influence Functional 

In this subsection, we make connection with the treatment of QBM based on the influence 
functional jHj with a phase space representation for the Wigner function|15j. First we con- 
sider the case without an environment. We define the transition element between the initial 
state |xo (fa) at £ = and the final state \x q) at time t to be 

K(x, q; t | x , q ; 0) = (x q\e~ iHt \x q ). (10) 

The Liouville equation for the density matrix is 

i^p(t) = [H, P (t)}, (11) 

where [, ] is the commutator. In the coordinate representation, the density matrix becomes 

p(x, x', q, q, t) = (x q\p{t)\x q) (12) 

with the collective notation for bath variables q = {q n }- The time evolution of the density 
matrix is given by 

,,,', t ) = f dvd&vWCfrat | x , 5ol 0)p(x , x;, ft, 0)K"(x l ,q'-, t | xi,,; ; 0),(13) 

In the present paper, we assume that the characteristic time scale for the bath is much 
shorter than the system. Under this condition, we may integrate out the bath harmonic oscil- 
lator variables to obtain an equation for the reduced density matrix p r (x, x') = j dqp(x, x', q, q, t). 
For a factorized initial condition between the system and the bath, which is assumed to be 
initially in thermal equilibrium, 

p(x , x' , g , q'o, 0) = ps(x , x' , 0) ® p B (q , q' , 0), (14) 
we can express the time evolution for the reduced density matrix in an integral form, 

p r (x,x',t) = J dxodx' J r (x, x'\ t | xo, x ; 0)ps(xo, x' , 0), (15) 

where its time evolution operator is given by 

J r (x, x'\ t | x , x' q ] 0) = J dqdq dq' Q K(x, q; t \ x , q ; 0)p B (g , q' , Q)K*(x', q; t | x' , q' Q ; 0) (16) 

For a harmonic oscillator bath, we have the exact expression 

J r (x, x'; t | x , x' ; 0) = / VxVx'e iS[x ' x ' ] . 

J(XQX' ) 

(17) 
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The total action iS[x,x'] consists of several contributions: 

<S[x,x'] = Ss[x,x'] + ASc[x,x'] + Se[x,x'] + Sip[x,x'], (18) 
where the sum of the actions for the system S plus its counter action is given by 

(S s + AS s )[R 1 r}= f ds{M R(s)r(s) - M nlR(s)r(s)}, (19) 
Jo 

where R=(x + x')/2, r = x — x' and for notational convenience, we have assumed the bare 
mass Mq and bare frequency Qq take on the values M = M and MqQ 2 = MQ 2 + AM'-fA/n 
for v = 1 while M = M + AM-fA/n and M Q 2 = MQ 2 + 4M7A 3 / (3ir) for u = 3. The action 
for the external field is 

S E \R,r] = C dsr(s)E(s). (20) 
Jo 

The influence action <S/ir[x, x'\ accounts for the effect of the bath on 5* and is given by 

Sjf[R, r] — % ds ds'r(s)fj,(s — s')r(s') 
Jo Jo 

-2 f ds f ds'r(s)u(s- s')R(s'), (21) 
Jo Jo 



where 



1 f°° (3Tiuj 

v(t) = — du J b(uj) coth coscut, (22) 

n Jo 2 

1 f°° 

n(t) = / duJE^uj) smut (23) 

7T JO 



are the noise and dissipation kernels respectively. 

From Eqs. (J21j) the Euler-Lagrange equations for R and r are 

M RJt) + M Q n 2 R c (t) + 2 / dsfj,(t - s)RJs) = E(t), (24) 

Jo 

M r c {s) + M Q 2 r c (s) - 2 J* ds'^(s - s')r c (s') = 0. (25) 

These equations have nonlocal kernels which contain the information of the past history of 
the bath in the presence of the system variables. Because of this, these equations normally 
contain time derivatives higher than two. As a result, they admit unphysical solutions. 
These unphysical solutions are removed by an order reduction procedure, reducing them 
into second order differential equations with well-defined initial value problems. They can 
also be specified uniquely by imposing the initial and final conditions: Ro and Rt (tq and 

n)- 
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If we let the two independent solutions of the homogeneous part of Eq. (jUj) (Eq. (j2*5|)) 
be Ui(s)(vi(s)), i = 1,2, with boundary conditions «i(0) = l,ui(t) = 0, w 2 (0) = 0,u 2 (t) = 1 
(t>i(0) = l,vi(t) = 0, ^2(0) = 0, V2(t) = 1), the solutions of these uncoupled equations can 
be written as 

R c (s) = RoUi(s) + Rtu 2 (s) + e(s), 

r c {s) = r vi(s) +r t v 2 {s), (26) 

where e(t) = f dsg + (t — s)E(s)/M. The solutions V\ and v 2 satisfy the homogeneous part 
of the backward time equation (J25J) and are related to U\ and u 2 by V\(s) = u 2 (t — s) and 
v 2 (s) = ui(t — s). The function g+(s) (g_(s)) also satisfies the homogeneous part of Eq. (J23|) 
(Eq. (j2SJ)) with boundary conditions «7±(0) = 0,g±(0) = 1. The solutions for g± for Ohmic 
and supra-Ohmic cases are given in Appendix A of [TSj. From these solutions u^ 2 and V\, 2 
can be determined. 

Since the potentials in our model are harmonic, an exact evaluation of the path integral 
can be carried out. It is dominated by the classical solution given in (|2fi)l . From these 
classical solutions, we write the action S[x, x'] as 

S[R c ,r c ] = (Mu x (t)R Q + Mu 2 (t)Rt)r t 

- (Mu 1 (0)R + Mu 2 (0)R t )r 

+ i(a n (t)rl + (a 12 (t) + a 21 (t))r r t + a 22 (t)rfj 

+ ei(t)r + e 2 (t)r t . (27) 

Here (ei(t), e 2 (t)) = e T = J * ds(v i(s), v 2 (s))e(s) and 

1 ft ft 

a ki(t) = 77 / ds / dsv k (s)/i(s - s')vi(s'), (28) 



2 Jo Jo 

for (A;, Z = 1, 2) contains the effects of induced fluctuations from the bath on the system 
dynamics. 

Using the results above, J r in Eq. (J 17)) can be written in the compact form, 

J r {R t ,r t ;t I R ,r ;0) = N(t)e lC , (29) 
where C = R T ur + zr T ar + e T r, (a)y = a^-, R T = (R , R t ) and r T = (r , r t ) 





' Mil 


«12 






V W21 


«22 / 





-Ul(0) Mi(t) 

-u 2 (0) u 2 (t) 



(30) 



2.3 QBM in the phase space representation 

The Wigner function is related to the density matrix by 

WJR, P,t) = — [ dre- iPr p r (R + r/2, R - r/2, £). (31) 
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The Wigner distribution function obeys the evolution equation 



W r (R t , P t , t) = J dR dP K(R t , P t ; t | R , P ; 0)W r (R , P Q , 0), (32) 
where K(R, P; t \ R , P ; 0) is defined by 

K{R, P; t | Ro, P ; 0) = J drdr e'^^ J r (R, r; t \ R , r ; 0) . (33) 
The propagator K for the Wigner function is given by 

K(R, P; t | Ro, P ; 0) = ^ / drdr j(-rr+P r +C)/h 

= N w (t) exp [ - SX^^SX}] , 

where N w (t) = N(t)/2^J\a\ and |a| is the determinant of a. The vector SX = X — (X), with 

* = (p-J- (34) 



and 



<*>= <s = * =± } * . (35) 



(P) J \ C12 C 22 / V ^o + ei y m 2 i V |u| m 22 y \ P + ei 
Here S is a matrix characterizing the induced fluctuations from the environment: 
s = ^_ / an ai2«2i - an«22 \ 

U21 I a 12«21 - anM 2 2 OnM22 - 2ai 2 M2lM22 + 022^1 J ' 



(36) 



At long times, fluctuations of the system are goverened by these terms as En — > ((A/?) 2 ), 
S 2 2^((AP) 2 ),S 12 = S2i^0. 

It is seen that this solution for the density matrix obeys non-Markovian dynamics in 
that the solution at a given time depends on its past history. Owing to the time dependent 
nature of their coefficients, despite its simple appearance, these equations are not easy to 
solve without approximations. The commonly used Markovian approximations may miss 
the essential features in the description of quantum/classical correspondence: it tends to 
underestimate the loss of quantum coherence because the rapid initial increase of diffusion 
coefficients is crucial for decoherence at low temperature in the strong coupling case. It is 
simply not a good approximation for a harmonic oscillator model with a generic spectral 
density. 
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3 Effective spin-boson model from QBM 
3.1 Dynamical level reduction 

We first illustrate our scheme of dynamical level reduction based on the harmonic QBM, 
which can be viewed as an infinite-level system in a bosonic environment 

N B 

h£° ] + Hf° ] = n^a + V2Q(a + a*) ^ c n q n (37) 

n=l 

This can be viewed as a limit of the finite iV-level system: 

N B 

H (N) + H (N) = QS + S . + ^ + 5 ,+ ) g ^ (3g) 

n=l 



when A" — > oo. Here we have absorbed V2fi by defining c n = \/2Vtc n . 

At finite temperature T, only those modes up to A^ ~ ksT/hQ are occupied. Thus at 
low temperature T ~ HQ, the effective number of levels of harmonic QBM is significantly 
reduced. In particular, at fc^T < hfl, we expect that the system is effectively reduced to 
two-levels: 

N B 

H f + HP = QS+S^ + (£ 2 - + 5+) £ c n q n . (39) 

n=l 

The formal correspondence is achieved by replacing the harmonic oscillator annihilation/creation 
operator a,a^ by the two-level pseudo spin annihilation/creation (Pauli) operator S^S^. 
The spin-boson model can be obtained by rewriting the Pauli operators as 



, n e 

Hf + Hf ] = Q(S Z 2 +-) + 5» ~ c n<ln- (40) 



n=l 



3.2 Fock states from phase space representation 

The correspondence between the Fock state representation for the pseudo spin qubits and 
the phase space representation is given as follows. First we write the density matrix in terms 
of the phase space variable as 

p{t) = f—X Q (z,z,t)e—e-^\ (41) 

J 7T 

where 

XQ (z,z,t) = Tr[p(t)e iza e^} (42) 
is a characteristic function for the Q representation fTHj fT7] . In a Fock space representation, 

Pkt(t) = f—XQ(z,z,t)(k\e- iSa 'e- iza \l), (43) 

J 71 
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Xq(z,z) is related to the characteristic function for the Wigner representation Xw(z, z) by 

\z\ 2 

X Q (z,z,t) = e 2 xw(z,z,t). (44) 

These characteristic functions are Fourier components of the phase space distribution func- 
tions, thus 



XQ (z,z,t)= I d 2 aQ(a,a)e im e iza (45) 
Xw {z,z,t)= I ' d 2 aW{a,a)e im e iza . (46) 



The characteristic function xc)(z,z,t) for the harmonic QBM evolved from the initial 
ground state has the following Gaussian form: 

Xq(z, z, t) = exp — a{t)z 2 — a{t)z 2 — 2b(t)\z\ 2 + iaif(t)z + icef(t)z (47) 

The time dependent coefficients a = c + a, b = C + S + l/4, and ctf in the above have 
their origins in the classical trajectory C of a damped harmonic oscillator given in Eq. (|35|) . 
the induced fluctuations S from the bath given in Eq. (|3fij) . and the external field E. The 
relations of these components are given as follows: 



8c = 


r 2 - 


C 2 

C u + Q C 12 — + 2i(CnCi2 + C21C22) 




8C = 


C 2 n + 


r 2 

Q 2 r 2 , °21 , r 2 
11 U 12 + ^2 + 22 




Aa = 




^22 




4S = 




^22 


(4 



a f(t) = ^l (n + i£)g + (t-s)E(8), (49) 



and 



2flio \ eft, 

where g+ satisfies the homogeneous part of the equation of motion in (|24jh 

From Eq. (|43|) we can directly evaluate the density matrix in the Fock representation at 
arbitrary quantum number. For instance, for an initial ground state, p(0) = |0)(0|, in the 
absence of an external field, the ground state and the first excited state population can be 
written as 

Poof*) = TF> (50) 

2[b(t) 2 - \a(t)\ 2 ] l/2 



and 



1 6(0 , x 

Pn(t) = — 51 

2[b(t) 2 -\a(t)\ 2 } 1/2 4[b{ty -\a{t)\ 2 f 2 
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Let us introduce the Pauli spin representation for the two-level system: 



(°x(t)) = Pai{t) + Pio{t) 
(a y (t)} = ip 10 (t) - ip 01 (t) 

(c?z(t)) = pn(t) - Poo(t). 



(52) 



We can express them by the variables defined in Eg .([47 )1 - ([49 )1 for arbitrary two- level spin 
initial states as follows: 



M*)> 



l 



+ 



A[b(ty-\a(t)\ 



—^{[{a x (Q))C 22 - (a y (0))nC 12 ] [Re a(t) - b(t)} 



Im a(t) 



(53) 



Mt)) 



-i 



4[b(t) 2 -\a(t)\ 2 f 2 



(^(0))^-K(0))C n 



[Re a(t) + b(t)] 



and 



M*)> 



[((T X (0))C 22 - (a y (0))QC 12 } Im a(t)} 

m 



(54) 



4 [6(t) 2 — |a(t)| 2 ] 3/2 
[6(t) 2 -|a(t)|f /2 I 2 2 [6(t)2 _ [a(i)| 2 ] 3/2 / 



(55) 



The leakage at time t is given by C(t) = 1— min Tr Pp r (t) , where P is the projection oper- 
ator onto the computational subspace and the minimization is taken over initial conditions. 
In our case, P = J2n=o,i \ n )( n \- The source of the leakage in our model is the transition to 
higher modes. The leakage is typically estimated by perturbative methods. However, the 
exact temporal evolution of this function is highly nontrivial as we will see below. Note that 
from the form of our effective Hamiltonian in (|4U|) . the behavior of coherence and population 
between our model and some others in the literature (for example, in [S]) are interchanged. 
They are related to each other by a change of basis. We can obtain similar expressions in the 
presence of an external field. We will examine this case in Sec. 4.1. In the Markovian limit, 
if the limit exists, two-level spin states become coupled nontrivially and obey optical-Bloch 
type equations [TH]. 

3.3 Limitations of other approximations 
3.3.1 Born-Markov approximation 

In this approximation, the bath correlation is neglected. This may be obtained as a limit 
of high temperature or slow system evolution in the Ohmic bath or white noise bath. For a 
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generic bath spectral density, however, there is no such limit. In a supra-Ohmic bath, the 
diffusion constant, when time-averaged for a long time, vanishes owing to the ultra-short 
time correlation. In a sub-Ohmic bath, it diverges owing to the long-range correlations. 
Only Ohmic spectrum gives the finite constant diffusion term. 

For weak coupling, the off-resonant counter-rotating terms in the interaction Hamiltonian 
are often ignored by invoking the rot at ing-wave- approximation (RWA). Although the use of 
RWA significantly simplifies the analysis, the dynamics under this approximation cannot 
capture the fast dynamics at time scales less than the natural time scale of the system. 
Furthermore, the spectrum of the Hamiltonian under RWA is found to be unbounded from 
below [in] • These features suggest that the range of validity of RWA is restricted to the leading 
order in the coupling constant only, where the counter-rotating terms do not contribute. 
After neglecting the counter-rotating terms from the two-level spin-boson Hamiltonian in 
(J2IJ), we obtain 

N 

H s + H SB = nS+S 2 + (5 2 - + S+) £ c n q n ^H s + H RWA = ttS+S 2 + ]T c n (S 2 bl + S'+6 n )(56) 

n=l n=l 

where b n = (uj n q n +ip n ) / ^2uj n are bath annihilation operators. In the presence of an external 
field, under the RWA, the reduced density matrix for the Hamiltonian obeys an optical Bloch 
equation. This case is commonly described in quantum optics text books. 



3.3.2 Born-Markov-RWA in multi-level-system (MLS) 



For comparison, we make the same Born-Markov-RWA in our E2L-SBM. Since the naive 
high temperature limit of the master equation obtained from QBM violates positivity 26J, 
we start from the master equation in the Lindblad form[T2j. For a particle initially in the 
Fock state p(0) = \k)(k\, the Q distribution function at time t has the following form: 



Q(a, a) 



7T [1 + n B (l - e-T*)] 



cxp 



\a\ 



1 



k\ 



l + n B (l -e-T*) 
\a\ 2 e-^ 



(n B + l)(l-e-*) 
l+n B (l-e-T*) 



(n B + !)(!- e-T*){l + n B (l - e -*)} 



where n B = l/(e" — 1) is a Planck distribution factor. For p(0) = |1)(1|, from 



(57) 
,and 



[l + n B (l-e-T*)]' 



-it 



+ n B (l - e 



(58) 



and 



1 + e 



-it 



+ 



2e~~< t 



l + n B (l-e-T*) ( i + nB (i_ e -7t) [1 + rz B (l - e^*)] 



(59) 
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3.3.3 Born approximation 

It is known that any master equation can be written in a time-convolutionless form|27j. 
However, the exact master equation in this form is still difficult to deal with. Most approaches 
based on the master equation invoke Born-approximation, then we obtain the tractable form, 
which can be solved exactly for simple systems or numerically for othersj^U]. The master 
equation under weak-coupling approximation may be suitable for describing the short time 
dynamics but tends to predict incorrect behavior for long times |HE]- Our nonperturbative 
results free from the weak coupling approximation is applicable to arbitrary time scales. 

4 Results and Discussions 
4.1 Results 

In Fig. 1, the populations and the leakage at T = 50[mK], Q = l[GHz], 7 = 0.1[GHz], 
and A = 100[Gi7^] are shown. The initial state is assumed to be the first excited state. At 
this temperature, the exact and the Markovian results agree at an intermediate time scale 
(around t = 10[ns]) but disagree at initial times. The slow oscillations in Fig. lb of the exact 
curve are from effects due to counter-rotating terms, while the fast oscillations are due to 
the frequency cutoff. The large leakage indicates that at this temperature, k^T > HQ, 2LS 
description is not a good picture. In Fig. 2, T = case is shown. There is a drastic difference 
in the entire time range shown in the figure. The exact result follows the quick decay at early 
times up to t ~ 10 — 50 [ps]. Late time decay rate asymptotically approaches the value given 
by the Markov approximation. The leakage is relatively large initially but negligibly small at 
late times. This indicates that only the lowest two levels are essentially populated except for 
the initial times, t < 5[ns]. The initial rapid decay of population is originated in the large 
initial leakage due to the transition to noncomputational subspace. This may be related 
to the initial large increment of the diffusion constant in the exact master equation at low 
temperature jS]- The total decay slows down as the leakage is suppressed at an intermediate 
time scale. In Fig. 3, the result for a supra-Ohmic environment is plotted. Compared to 
the Ohmic case, the initial decay of the excited state population is much more drastic but it 
appears to saturate at late times. Thus if the initial decay is strong enough, the coherence 
in the system can be totally washed out at an early stage, a serious concern for the quantum 
devices. On the other hand, if it is small, the system can remain coherent for a long time. 
Note that our result disagrees markedly with the Markovian prediction over the entire time 
range. 

In Fig. 4, the Rabi oscillations in the presence of an external sinusoidal pulse at the 
resonant frequency are plotted. The most notable difference between the exact results from 
the Markovian results is that the exact results show the low onset and low visibility for all 
times. The difference is more evident for the supra-Ohmic case. Our figures also suggest that 
it is not easy to determine the characteristics of the environment only from the experimental 
Rabi oscillation data without the precise knowledge of the dissipation. The large increase 
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of leakage is due to the resonant transition to high level states. Though increasing anhar- 
monicity in the potential will suppress these transitions to some extent, the initial rapid 
increase of leakage is unavoidable due to energy-time uncertainty relation. In the presence 
of tunneling with a biased potential, due to the existence of resonant transitions to the con- 
tinuum modes, we expect the result will be qualitatively similar to our case. In this case, 
the leakage in our model can be interpreted as the effect of tunneling, or more appropriately, 
environment-induced hopping to other metastable states. 
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4.2 Discussion 



We saw that at low enough temperatures, many conventional approaches based on the Born- 
Markov approximation can significantly underestimate the environment-induced decoherence 
beyond the weak system-bath coupling. In this regime, the visibility in Rabi oscillations 
in the exact calculation tends to be lower than what is expected in the Markovian ap- 
proximation. Low visibility in Rabi oscillations is commonly observed in superconducting 
qubits[2~Tl |22 I2H1 121] ■ The bath time scale is also important in causing the initial rapid 
decoherence and leakage; this is completely neglected in analysis based on the Born-Markov- 
RWA. This initial effect can manifest itself as an onset value of Rabi-oscillations. In many 
practical implementations of qubits, the temperature of the environment compared to the 
bath is small, ksT << hA, thus we are in the low temperature regime. 

The E2L-SBM approach gives a precise evaluation of the leakage due to the system's 
interaction with the environment and the external control field. For temperatures higher 
than the characteristic energy of the oscillator, the large leakage makes the qubit based on the 
choice of the lowest two levels ill-defined. During gate operations, this can become a serious 
problem and remedies for stabilizing the system such as using external pulse control may be 
necessary. Our result shows that the time scale associated with leakage is characterized by 
the dynamical time scale of the system and the bath. 

In realistic macro- or mesoscopic systems, the potential contains anharmonicity, which 
causes the deviation of the system dynamics from the harmonic motion. A measure of 
anharmonicity near ground states can be given by the ratio of energy level separation between 
the lowest levels Uqi = iO\ — u and the excited levels ujyi = uo 2 — When this ratio is 
small, cg>i 2 /cg>oi << ^01 , the initial short time evolution around a metastable state can be 
well- described by the linear dynamics. When the correction to the energy level due to 
anharmonicity in the potential becomes important, it is necessary to include such an effect 
in our scheme. Although the large anharmonicity also prevents the leakage in the long term, 
the initial large leakage we saw cannot be completely eliminated for the reason we mentioned 
before. When we apply our formalism to the metastable state, eventually the system state 
will leave the harmonic oscillator phase space into other metastable states via tunneling. 
The harmonic approximation of coherent dynamics is expected to be accurate at initial 
times when the time scales associated with these nonlinear effects are large compared to 
the decoherence time scale. In our example, the deviation from the Markovian prediction is 
evident in the very early stage of the system evolution up to t ~ 1 [ns] even for an intermediate 
temperature. For the realistic implementation of qubits, the underlying potential landscape 
leading to the discrete energy level is already known by design[2H 1221 El I2H EHl ED] - Our 
approach based on E2L-SBM is suitable in this situation and will give a more precise estimate 
of the open system dynamics than the one based on the conventional 2LS approximation. In 
particular, our results are directly relevant to the superconducting phase qubit models [2U 
1221 12H] ■ in the superconducting qubits, the major source of decoherence is the noise induced 
by the interaction with the current or charge sources mainly during the qubit manipulations. 
We have not considered other possible sources of decoherence such as the coupling to defects 
or nuclear and magnetic spins. Multilevel structure in the superconducting flux qubits was 
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studied in [3U| by Born-Markov approximation without control fields. 

For the system-environment coupling we considered in Eq.Q, the Fock state is not 
an eigenstate of the interaction Hamiltonian and is subjected to a complex decay even 
under the Born-Markov approximation as shown in Sec. 3.3.2. Previous study in the high 
temperature limit indicates the pointer state under this system-environment coupling is a 
coherent state [7j. Our calculation based on the exact solution for QBM indicates that, 
beyond the weak coupling regime, the environment-induced effect has a crucial impact on 
the system dynamics at an early stage. A factorized initial condition is used to derive our 
main results in accord with the initialization scheme used commonly in quantum information 
processing [HI]. From the decoherence study of QBM in the presence of the initial system- 
bath correlation due to the preparation effect |U E2] , we expect that our results are robust 
and should hold for a more general class of initial conditions. 
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Figure 1: (a) Plot of the time-evolution of the population of the ground and the first excited 
state population, and the leakage (the dot-dashed curve) at T = 50[mfC] with Q = l[GHz], 
7 = 0.1[GHz], A = 100[GHz]. The thick solid (dashed) curve is the exact (Markovian) 
result for excited states while the thin solid (dashed) curve is the exact (Markovian) result 
for ground states. Panel (b) is in the logarithmic time scale. 
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Figure 2: The ground and the first excited state population, and the leakage (the dot-dashed 
curve) at T = 0. Q = l[GHz], 7 = 0.1[GHz], A = 10O[GHz]. The thick solid (dashed) 
curve is the exact (Markovian) result for excited states while the thin solid (dotted) curve is 
the exact (Markovian) result for ground states. Panel (b) is in the logarithmic time scale. 
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Figure 3: The excited state populations and the leakage are plotted for supra-Ohmic envi- 
ronment. T = 10[miT] in panel (a) and T = 50[miT] in panel (b). The thick (thin) solid 
curves is an exact result for supra-Ohmic (Ohmic) case, while the dashed line is from the 
Born-Markov approximation. Q = 1.5[GHz], 7 = 0.1, A = 100[GHz\. 
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Figure 4: The excited state populations in the presence and absence of external field at T = 
10 [mK] are plotted for Ohmic environment in panel (a) and for supra-Ohmic environment 
in panel (b). The applied external is sinusoidal in the form: E = Eocos(Qt) with Eq = 1.0. 
The thick (thin) curves are exact results in the presence (absence) of external field, while 
the dashed curves are from the Born-Markov apporoximation. Q = 1.5 [Gifz], 7 = 0.1, A = 
100[GHz}. 
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